Mon. Not. R. Astron. Soc. Q00,[T|fl6l(2009) Printed 22 May 2009 (MN KTeX style file v2.2) 



HiZELS: a high redshift survey of Ha emitters. II: the nature of 
star-forming galaxies at z=0.84 ^ 

0^ D. SobralH, P. N. Best\ J. E. Geach^, Ian Smail^, J. Kurk^, M. Cirasuolo\ M. Casali^ 
§ R. J. Ivison^ ^ K. Coppin^ & G. B. Dalton^'^ 

^SUPA, Institute for Astronomy, Royal Observatory of Edinburgh, Blackford Hill, Edinburgh, EH9 3HJ, UK 

^Institute of Computational Cosmology, Durham University, South Road, Durham, DHl 3LE, UK 

^ Max-Planck-Institut fur Astronomie, Konigstuhl, 17 D-69117, Heidelberg, Germany 

^European Southern Observatory, Karl-Schwarzschild-Strasse 2, D-85738 Garching, Germany 

^Astronomy Technology Centre, Royal Observatory of Edinburgh, Blackford Hill, Edinburgh, EH9 3HJ, UK 

^Astrophysics, Department of Physics, Keble Road, Oxford, 0X1 3RH, UK 

^ Space Science and Technology, Rutherford Appleton Laboratory, HSIC, Didcot, 0X11 OQX, UK 



Accepted 2009 May 22. Received 2009 May 11; in original form 2009 January 25 



Oh 



ABSTRACT 

New results from a large survey of Ha emission-line galaxies at z = 0.84 using the Wide 
Field CAMera on UKIRT and a custom narrow-band filter in the J band are presented as 
O part of the HiZELS survey. The deep narrow-band images reach an effective flux limit of 

^ Fhq! ~ 10~^^ergs~^ cm~^ in a co-moving volume of 1.8 x 10^ Mpc^, resulting in the largest 

^ and deepest survey of its kind ever done at 2; ~ 1. There are 1517 potential line emitters 

■ detected across ~ 1.4 square degrees (in two fields: COSMOS and UKIDSS UDS), of which 

743 are selected as Ha emitters, based on their photometric and spectroscopic redshifts. These 
are then used to calculate the Ha luminosity function, which is well-fitted by a Schechter 
> function with L* = io42-26±o.05 ergs-\ (^*=10-^-^2±0-^o Mpc"^ and a = -1.65 ± 0.15, 

and are used to estimate the volume average star formation rate at z = 0.845, psfr' 
0.15 ± 0.01 M0 yr~^ Mpc~^ (corrected for 15% AGN contamination and integrated down 
to 2.5 M0yr~^). These results robustly confirm a strong evolution of psFR from the present 
day out to z ~ 1 and then flattening to z ~ 2 using a single star-formation indicator: Ha 
^-H luminosity. Out to z ~ 1, both the characteristic luminosity and space density of the Ha emit- 

ters increase significantly; at higher redshifts, L* continues to increase, but decreases. The 
z = 0.84 Ha emitters are mostly disk galaxies (82 ± 3%), while 28 ± 4% of the sample show 
signs of merger activity; mergers account for ~ 20% of the total integrated Psfr at this red- 
^ shift. Irregulars and mergers dominate the Ha luminosity function above L*, while disks are 

dominant at fainter luminosities. These results demonstrate that it is the evolution of "normal" 
rN disk galaxies that drives the strong increase in the star formation rate density from the current 

^ epoch to 2: ~ 1, although the continued strong evolution of L* beyond z = 1 suggests an 

increasing importance of merger activity at higher redshifts. 

Key words: galaxies: high-redshift, galaxies: luminosity function, cosmology: observations, 
galaxies: evolution. 
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1 INTRODUCTION 

Understanding the basic features of galaxy formation and evolu- 
tion requires unveiling the volume- averaged star formation rate as a 
function of epoch, psFR, its distribution function within the galaxy 
population, and its variation with environment. At the present, sur- 
veys of the star-formation rate density suggest a rise of (1 + z)^ 



* Based on observations obtained with the Wide Field CAMera (WFCAM) 
on the United Kingdom Infrared Telescope (UKIRT) 
t E-mail: drss@roe.ac.uk 



out to at least z ^ 1 (e.g.pUy et al.|1996[ [Connolly et al.|1997[ 
Hogg et al.|1998l|Haarsma et al.|2000[|Hopkins|2004[ [Hopkins &| 
Beacom||2006| l indicating that most of the stars in galaxies today 
formed at 2: > 1. 

In the local Universe, studies have demonstrated that star for- 
mation is strongly dependent on the environment. While clusters 
of galaxies seem to be primarily populated by passively-evolving 
galaxies, star-forming galaxies are mainly found in less dense envi- 
ronments (e.g. [Lewis et al.|2602 ). Star-forming galaxies have also 
been found to have lower masses than passive galaxies ( [Cowie etaL] 
[1996| ). How do these environmental and mass dependencies change 
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with cosmic time? When did they start to be noticeable, and how do 
they affect the evolution of galaxies, clusters and the star formation 
rate density of the Universe as a whole? How much of the evo- 
lution of the cosmic star formation rate density is associated with 
the evolution of star-forming galaxies and how much is driven by 
galaxy merger activity? In order to properly answer such questions 
it is mandatory to conduct observational surveys at high redshift, 
which can then be used to test theoretical models of galaxy evolu- 
tion. These can be performed with a large variety of techniques and 
instrumentation. 

There are many different star formation indicators which have 
been widely used both to select and study star-forming galaxies, 
such as the Ultra- Violet (UV) light emitted by young stars, the en- 
ergy that is absorbed and then re-emitted in the form of emission 
lines (such as Ha or [Oil] 3727), or in the sub-mm and infra-red, 
or even radio emission from supernovae. Ideally, the use of dif- 
ferent star formation tracers would provide consistent answers to 
most of the questions that remain unanswered. Unfortunately, stud- 
ies with different star formation tracers suffer different biases and 
selection effects, which contribute to considerable discrepancies. 
These problems are also amplified by the effects of cosmic variance 
in most of the current samples. Another issue is the difficulty in cor- 
recting for extinction, especially for UV and optical wavelengths, 
which can lead to large systematic uncertainties in the star forma- 
tion densities derived from measurements in these wavebands. 

No single indicator provides a perfect view of the evolution of 
the star-formation rate density ( K ewley et al._200 4 ), but using dif- 
ferent tracers is not the solution either. Ideally, one needs to choose 
a star-formation indicator that can be applied from the present to 
high redshift, is relatively immune to dust extinction, and which has 
sufficient sensitivity to ensure that the derived star-formation rate 
distribution does not require large extrapolations for faint sources 
below the sensitivity limit. The Ha luminosity is arguably the best 
candidate to achieve these goals and large-area surveys with a sen- 
sitivity of a few solar masses per year can now be done with current 
instrumentation out to 2; ^ 2. This contrasts with the equivalent 
star formation rate limits for dust-independent tracers such as ra- 
dio, far-infrared or sub-millimeter observations which vary from 
- 10-100 Mo yr"^ at the same redshifts for large area surveys. 

Large Ha surveys have been carried out since the 1970's (e.g. 
Cohen|1 976, Kennicutt & Kent 1983, Romanishin 1990, Gavazzi 
etal.|1991 , Galleg o et al.|1995HYan et al.|1999„Tresse et al.,2002t 



Doherty et al. 200^. With the rise of large-format imaging cam- 



eras in the optical and in the near-infrared, narrow-band filters can 
now be used to undertake deep surveys for emission-line objects 
in large volumes which cover the great majority of the cosmic his- 
tory, out to 2; 2.5 using Ha (e.g. |Villar et al.|2008[[Geach et al.] 
|2008| ). The surveys identify sources on the basis of the strength of 
their emission line and thus represent roughly a star-formation rate 
selected sample. Hence, by using a set of narrow-band filters, it is 
possible to apply a single technique to target Ha emitters across 
a wide range of redshifts, gathering representative samples at each 
epoch with a uniform selection. Those samples are powerful tools 
for tracing the evolution in the star-formation rate density across 
the expected peak of star-formation in the history of the Universe, 
which is one of the main aims of the High- 2; Emission Line Survey 
(HiZELS). 

HiZELS is a panoramic extragalactic survey using the wide 
area coverage of the Wide Field CAMera (WFCAM, Casali et al. 
2007) instrument on the 3.8-m UK InfraRed Telescope (UKIRT). 
The survey utilizes a set of existing and custom-made narrow-band 
filters in the J, H and K bands to detect emission line galaxies 



over 5 square degrees of extragalactic sky. The H2S(1) narrow- 
band filter is being used to target Ha emitting galaxies at 2; = 2.23, 
which are predominantly star-forming systems (Geach et al. 2008). 
In addition, narrow-band filters in the J and H bands (hereafter 
NB J and NB//) have been custom-designed to target the [Oil] 3727 
and [Olll] 5007 emission lines in galaxies at the same redshift as the 
H2S(1) Ha survey. Together, the three sets of filters help to detail 
the properties of the line emitters at ^ = 2.23, while the NB j and 
NBi/ filters deliver identically- selected Ha samples at 2; = 0.84 
and 1.47 respectively. 

This paper presents deep narrow-band imaging using the NB j 
filter at A = 1.211|j,m, as part of HiZELS, over ^ 0.7 deg^ 
of the in the SX DF Subaru-XMM-UKIDSS Ultra Deep Survey 
( [Lawrence et~ari[2007| ) field (UDS), and - 0.8 deg^ in the Cos- 
mological Evolution Survey ( [Scoville et al.||200"7{ [Koekemoer et| 
|al.|2007| ) field (COSMOS). This corresponds to an area coverage 
which is ~ 8 times larger and twice the depth of the survey by| |Vil-| 
|lar et al.| ( |2008| ), the best previous emission-line survey ai z ^ 1. 

The paper is organised in the following way. §2 outlines the 
details of the observations, and describes the data reduction, photo- 
metric calibration, source extraction and survey limits. §3 presents 
the narrowband selection criteria, the final sample and the photo- 
metric and spectroscopic redshift analysis. Results are presented in 
§4: the Ha luminosity function, the star-formation rate density, the 
morphologies of the Ha emitters and their relation and contribution 
to the Ha luminosity function. Finally, §5 outlines the conclusions. 
An Ho — 70kms~^ Mpc~\ Qm — 0.3 and Qa - 0.7 cosmology 
is used and all magnitudes are in the Vega system, except if noted 
otherwise. 



2 OBSERVATIONS AND DATA REDUCTION 

Both fields were observed with WFCAM on UKIRT using a set 
of custom narrow-band J filters (A = 1.211iJ,m, = 0.014iJ,m), 
during 2007 January 13-17 and February 13,14 & 16 for the COS- 
MOS field and October 19-24 for UKIDSS UDS. WFCAM's stan- 
dard "paw-print" configuration of four 2048 x 2048 0.4^^ pixel" ^ 
detectors offset by ~ 20^ can be macrostepped four t imes to cover a 
contiguous region of ^ 55^ x 55' ( [Casali et ah 2007| . For each field 
(COSMOS and UDS) 4 paw-prints, or ^ 0.8 deg^ were mapped 
with narrowband exposures of ~ 2 Iks pixel" ^. The seeing varied 
between 0.8-1.0^' during the observing nights, and conditions were 
photometric. The NDR (Non Destructive Read) mode was used for 
all narrow-band observations to minimise the effects of cosmic rays 
in long exposures. The observations were obtained following a 14- 
point jitter pattern for UDS (to match the UKIDSS observing strat- 
egy) and a 9-point jitter sequence for COSMOS. Broad-band J- 
band observations were also obtained in the COSMOS field, but 
were not necessary for the UDS due to the availability of UKIDSS 
data. A summary of the observations is given in Table [T] 



2.1 Data reduction 

In order to accomplish the survey goals, a dedicated pipeline has 
been developed: PfHiZELS. This was done following the method 
presented in |Geach et al.| (|2008 ), improving several steps and pro- 
ducing a user-friendly, completely automated pipeline. It consists 
of a set of PYTHON scripts which implement a complete reduc- 
tion of WFCAM data, taking advantage of several publicly avail- 
able packages, as detailed below. Small samples of data were also 
reduced by hand using IRAF (including DIMSUM and XDIMSUM), 
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Table 1. Observation log for the NB j and J observations of the COSMOS and UKIDSS UDS fields. 
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iv./v. 
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Int. time 
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Dates 


'^lim 










(Ks; 


K ) 


i^ZUU / ) 


per; 


COSMOS 1 


NBj 


100000 


+021030 


19.7 


1.0 


14-16 Jan 


21.7 


COSMOS 2 


NBj 


100052 


+021030 


21.6 


1.0 


13, 14 Jan 


21.6 


COSMOS 3 


NBj 


100000 


+0223 44 


19.0 


0.9 


15-17 Jan 


21.7 


COSMOS 4 


NBj 


100053 


+0223 44 


17.2 


0.9 


15, 17 Jan, 13-16 Feb 


21.6 


COSMOS 1 


J 


100000 


+021030 


5.7 


0.9 


14-16 Jan 


22.8 


COSMOS 2 


J 


100052 


+021030 


6.9 


1.0 


13, 14 Jan 


22.8 


COSMOS 3 


J 


100000 


+0223 44 


5.7 


0.9 


15-17 Jan 


22.8 


COSMOS 4 


J 


100053 


+0223 44 


5.1 


0.9 


15, 17 Jan, 13-16 Feb 


22.6 


UDSNE 


NBj 


0218 29 


-045220 


21.2 


1.0 


18, 20 Oct 


21.6 


UDSNW 


NBj 


0217 36 


-0452 20 


22.6 


1.0 


19, 21 Oct 


21.7 


UDS SE 


NBj 


0218 29 


-05 05 53 


20.0 


1.0 


19 Oct 


21.5 


UDS SW 


NBj 


0217 38 


-05 05 34 


22.5 


1.0 


20, 21 Oct 


21.6 



which allowed a comparison with the automatic data reduction and 
showed that the results are comparable. 

Dark frames were median combined for each observing night 
and the individual science frames were then dark subtracted. Next, 
frames were median combined in each jitter sequence of 14 (UDS) 
or 9 (COSMOS) without any offset, to produce a rough first pass 
flat. The latter was then used to produce a badpixel mask for each 
chip by flagging pixels which deviate by more than 3a from the 
median value. Each frame was then flattened, and an individual 
badpixel mask (for each frame) was produced. 

In order to produce better flats, S Extractor ( [Bertin &] 
|Arnouts|1996| ) was used to identify objects on the flattened images, 
and source masks were created and used to produce second-pass 
flat fields. These were then used to create the final flattened frames. 

A world coordinate system was fitted to each frame by query- 
ing the USNO A2.0 catalog. This fitted, on average, ^ 100 ob- 
jects per frame and rotated images appropriately for each WFCAM 
chip. Frames were then de -jittered and co-added with SWarp 
(Bertin 1998), which performs a background mesh-based sky sub- 
traction optimized for narrow-band data. It should be pointed out 
that WFCAM frames suffer from significant cross-talk artifacts, 
which manifest themselves as toroidal features at regular (every 
128) pixel intervals from sources, in the read-out direction. Fur- 
thermore, as these are linked to a "physical" location, they can only 
be removed from the source catalog (see §2.3| ). 



2.2 Photometry calibration 

Narrow- and broad-band images were photometrically calibrated 
(independently) by matching ~ 70 stars with mj = 11-16 per 
frame from the 2MASS All-Sky catalog of Point Sources ( |Cutri| 
let al.|2003j which are unsaturated in narrow-band frames. Based 
on the photometric calibration, the absolute calibration is expected 
to be accurate to < 3% when compared to 2MASS. We have also 
double-checked the offset between zero-points by comparing them 
to the expected values, roughly given by —2.5 log (SAjv^^^x^t^^ ' 
where (5 A is the filter width and 6t the exposure time per frame. For 
convenience, after this step, narrow and broad-band images were 
normalised to give them the same zero-point. 



2.3 Source Extraction and Survey Limits 

The survey is made up of a mosaic of eight WFCAM pointings (4 
each for UDS and COSMOS), i.e., 8 x 4 x 13.7' x 13.7' tiles. 
The UKIDSS UDS J-band image does not overlap entirely with 
the full narrow-band image, leading to a total overlapping area of 
0.70 deg^ while for COSMOS the overlap is 0.78 deg^ 

Sources were extracted using S Extractor ( [Bertin &| 
[Arnout s 1996 ). Optimal parameters were found by running a large 
number of different extractions which converged to a set of parame- 
ters that allowed the extraction of all obvious sources down to the 3- 
a limit in each frame and minimized the extraction of noise/artifact 
features. The extraction included an optimized sky subtraction, and 
fixed photometry apertures of 3'' (diameter) were used. Several 
tests were done using S Extractor in dual mode (using frames 
in one band to detect the sources and measuring on the other band) 
and single mode. The first option was used for COSMOS (where 
data were taken in both bands and so the frames were extremely 
well registered), while for UDS (where UKIDSS data was used for 
the J band) the extraction was done on both bands independently, 
followed by a match using a simple 2.0" criteria. In 11 cases there 
was more than one match; for these, a careful local astrometry so- 
lution was calculated based on nearby sources clearly identified in 
both bands and the correct source was then clear in all cases (the 
correct match was always the original closest). Furthermore, tests 
were run on COSMOS frames to show that both options produce 
comparable samples, and so this is not likely to produce any sig- 
nificant difference between the extraction for UDS and COSMOS. 
Narrowband sources with no clear J band detection were retained 
but there are few of these as the J coverage is significantly deeper 
than the NB j . As the observations have slightly different FHWM, 
different total exposure times, and because each WFCAM chip has 
slightly different properties, objects were extracted down to each 
chip's limit, which was then confirmed using Monte Carlo simula- 
tions (see §4.1.2| ). 

In order to clean spurious sources from the catalogue (es- 
sential to remove cross-talk artifacts), the frames were visually 
inspected showing that sources brighter than ^ 12 mag are sur- 
rounded by large numbers of artifacts detected within "bright ha- 
los", as well as cross-talk, while fainter sources (up to 16 mag) 
show only cross-talk features. Sources fainter than 16 mag do not 
produce any detectable cross-talk at the typical depth of the ob- 
servations. Cross-talk sources and detections in the halo regions 
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Figure 1. Two narrow-band excess objects selected from the NB j imag- 
ing with different colour significances (upper panels), and a strong isolated 
emitter (lower panels). Images are ~ 20^^ x 20^^ 



were removed from the catalogue separately for each frame, which 
greatly simplifies their identification (Geach et al. 2008). 

Catalogue sources detected in regions with less than 85% of 
the total integration time were removed; which still assures a com- 
plete overlap between frames. When a source was catalogued in 
more than one final image, the catalogue entry with a higher expo- 
sure time was selected. 

The average 3-cr depth of the entire set of NBj frames is 
21.6 mag, with J depths being J = 23.4 mag (UKIDSS UDS DR3) 
and J = 22.7 mag (COSMOS). The narrow-band imaging detects 
a total of 21773 objects in COSMOS across 0.76 deg^ and 15449 in 
UDS across 0.68 deg^. These areas include the removal of regions 
in which cross-talk and other artifacts caused by bright objects are 
located. 



3 SELECTION 

3.1 Narrowband excess selection 

Emission line systems are initially selected according to the signifi- 
cance of their (J— NB j) colour, as they will have (J— NB j) > 0. 
However, true emitters need to be distinguished from those with 
positive colours due to scatter in the magnitude measurements and 
this is done by quantifying the significance of the narrowband ex- 
cess. The parameter H quantifies the excess compared to the ran- 
dom scatter expected for a source with zero colour (Bunk er et al.| 
tl995). In other words, a source can be considered to have a genuine 
narrow-band excess if: 



CnB r — Cj > 



(1) 



where cnbj and cj are the counts for the NBj and J bands, re- 
spectively, and 6 is the combined photometric error: 



(2) 



The approach is similar to a standard signal-to-noise selection. 
Colour and E significances are related by 



J - NBj = -2.51og(l - E510-° '*(^^-'~'''^>), 



(3) 



where ZP is the zeropoint of both frames, assuming that the sky 
variation is the dominant feature contributing to the errors in pho- 
tometry. This has been tested by measuring counts in randomly 



placed apertures (~ 10000 for each frame), confirming the hypoth- 
esis. 

As the NBj filter does not fall at the centre of the J-band, 
objects with redder colours will tend to have a negative ( J— NB j) 
colour, while bluer sources will have (J— NBj) > 0. This will 
affect the selection of emission line objects, but can be corrected 
for by considering the broad-band colours of each source. To do 
this, (J-NBj) was plotted as a function of {z - J) (COSMOS) 
or {J — K) (UDS) colour and a linear fit was derived, determin- 
ing trends with a slope of (J-NBj) - 0.06(J - K) for UDS 
and ( J- NBj) - 0.11(z - J) for COSMOS. This was then used 
to correct the NBj magnitudes and thus the resultant (J— NBj) 
colours. Sources with no z or K data - either because they are 
too bright or too faint - were corrected assuming a median colour 
(based on sources at similar magnitudes). 

The flux of an emission line, Fune, and the equivalent width, 
EWiine, can then be expressed as: 



Sne 



fj 



{AXnbj/AXj 



and 



EWii 



■ Snbj 



fj 



fj-fNBj{AXNBj/AXj) 



(4) 



(5) 



where AXnbj and AAj are the widths of the broad and narrow- 
band filters, and fNBj and fj are the flux densities measured for 
each band. EWune is simply the ratio of the line flux and continuum 
flux density. 

The selection of emission-line candidates is done imposing 
two conditions. Firstly, the NB j sources are considered candidate 
emitters if they present a colour excess significance of S > 2.5. 
Secondly, only NBj detections with EWune > 50A (correspond- 
ing to (J— NBj) > 0.3) are selected. This is done to avoid in- 
cluding bright foreground objects with a large significance and a 
steep continuum across the J band (e.g. |van der Werf et al.|2000"| , 
and was chosen to reflect the general scatter around the zero colour 
a bright magnitudes. A comparison with spectroscopic data (see 
§3.5| ) indicates that these criteria maximise the completeness of the 
sample without introducing a significant number of false emitters. 
Figure [2] shows the colour-magnitude diagrams with the selection 
criteria for COSMOS and UDS, respectively. Figure[T]presents two 
examples of emitters. 



3.2 Complete sample 

Narrow-band detections below the 3-cr threshold were not con- 
sidered. When there was no > 3a J-band detection, a 3-cr up- 
per limit for the (J— NBj) colour was computed, as shown in 
Figure [2l The average 3-cr line flux limit over both fields is 
8 X 10~ ^ergs~^cm~^. The complete sample has 1517 excess 
sources out of all 37222 NBj detections in the entire area, with 
824 being detected in COSMOS (corresponding to 1084 emitters 
per deg^) and 693 in UDS (1020 emitters per deg^). The poten- 
tial emitters were all visually inspected in both J-band (when de- 
tected) and NBj, and 48 were removed from the sample as they 
were flagged as spurious. The majority of these (34) correspond to 
artifacts caused by bright stars that are on the edges of two or more 
frames simultaneously. The remaining 14 sources removed were 
low S/N detections in noisy regions of the NB j image. 
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Figure 2. Colour-magnitude plots for COSMOS (left panel) and UDS (right panel) showing all > 3a detections in the NB j image. The curves represent 
E significances of 5, 3, 2.5 and 2, respectively, as defined in the text (for an average depth). The dashed line represents an equivalent width cut of 50A. All 
selected narrow-band emitters are plotted in black, while candidate Ho; emitters (selected using photometric redshifts from Mobasher et al. 2007 for COSMOS 
and Cirasuolo et al. 2008 for UDS) are plotted in red. Emitters with no clear J-band detection are assigned a 3-cr upper limit. 



3.3 Photometric redshift analysis 

In order to select Ha emitters ai z = 0.84, one needs to separate 
them from other line emitters at different redshifts. The detection 
of H;^ and [Olll] 5007 emitters at z ^ 1.4-1.5 is expected, to- 
gether with [On] 3727 at 2; = 2.23 and other emission lines. Multi- 
wavelength data, photometric and spectroscopic redshifts available 
for both COSMOS ([Mobasher et al.pOOTj [Lilly et aL||2007 ) and 



UDS ( [Cirasuolo et al.|20 08|) are therefore used to distinguish be- 
tween different emission lines and also to evaluate how robust the 
emitter selection criteria are. In order to do this, photometric red- 
shifts for the emission line candidates were taken from the photo- 
metric redshift catalogues for COSMOS and UDS ( [Mobasher et al.[ 
|2007l [Cirasuolo et al. 2008 ). For UDS, the overlap with the photo- 
metric redshift catalogues reduces the area coverage to 0.54 deg^ 
due to the required overlap with Subaru optical data and further 
masking around bright stars: 133 sources are "lost" as they are in 
the excluded area. In addition, for 14 sources photometric redshifts 
are not available in the catalogues due to i/Ks magnitude limits (3 
in COSMOS, 1 1 in UDS): while these appear to be real sources (al- 
though there is a chance they are spurious), they are likely to be at 
redshifts higher than z = 0.84 and so are excluded from this analy- 
sis of Ha emitters. The photometric redshift matches result in a fi- 
nal sample of 1370 excess sources (821 in COSMOS, 549 in UDS) 
meeting all of the selection criteria across a total of 1.30 deg^. 

The photometric redshifts in COSMOS are particularly good 
for selecting Ha emitters at 2; = 0.84 as the available data are 
ideal to probe the 4000A break at that redshift. The deeper near- 
infrared data used in the UDS photometric redshifts provide better 
accuracy for higher redshift sources, and are accurate enough to 
probe the proportion of 2; ^ 1.4-1.5 [Olll] 5007 and H/3 emitters. 
Photometric redshifts for COSMOS present a{Az) = 0.03, where 



Az = {zpho 



). The fraction of outliers, defined 



as sources with Az > 3cr(Az), is lower than 3%, based on results 
from ^-COSMOS (see ^33). For UDS, the photometric redshifts 
have a{Az) = 0.04, with 2% of outliers. 

Figure [3] shows the photometric redshift distribution for the 
selected narrow-band emitters in COSMOS and UDS, demonstrat- 
ing very good agreement between the two fields, despite the com- 



pletely different photometric analyses - done with different codes, 
bands and by completely independent teams. Both samples peak 
dX z ^ 0.85, implying that the majority of the narrow-band ex- 
cess sources are indeed Ha emitters. In addition to this, there is an- 
other peak 3i z ^ 1.4-1.5 in UDS and the same - but at a slightly 
higher photo-z - in COSMOS. This is interpreted as a significant 
population of [Olll] 5007 emitters at z ~ 1.4 and/or Hf3 emitters 
ai z ^ 1.5. However, the photometric redshifts are not accurate 
enough to clearly distinguish between these two populations, and 
thus the slight difference between UDS and COSMOS might be just 
caused by the use of different bands and methods giving different 
performance at z > 1. 

3.4 Selecting Ha emitters ai z — 0.84 

The sample of candidate Ha emitters at 2; = 0.84 is derived by con- 
sidering not only the best-fit photometric-redshift for each source, 
but also more extensive information contained within the photo- 
metric redshift probability distribution. In particular, the Ha candi- 
dates are defined to be those sources with Zmin < 0.845 < Zmax 
(where Zynin and Zmax are the 1-cr redshift limits of the principle 
peak in the photometric redshift probability distribution). For COS- 
MOS, the photometric redshift selection produces a sample of 477 
potential Ha emitters, while for UDS the same procedure yields 
270 sources (see §3.3[ for details). 

While it is true that this selection can potentially introduce 
some biases, a careful comparison of the selected Ha emitters 
and the remaining emitters shows no evidence for such an effect. 
The emission line equivalent width distributions of the selected 
Ha and non-Ha emitters are statistically indistinguishable. The 
sample of non- selected emitters does contain objects which have 
slightly fainter magnitudes (and correspondingly fainter emission 
lines) than the selected Ha emitters, but the colours of these fainter 
galaxies are clearly consistent with them simply being at higher 
redshift (e.g. [OIII], H/3 or [Oil] emitters): they are typically de- 
tected at a ~20-30 a (or higher) level in the optical images, and 
thus it should not be any inaccuracy in the photometric redshifts 
which affects their selection. One can also see that the photomet- 
ric redshift distribution in Figure [4] (left panel) is clearly different 
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Figure 3. The left panel shows the photometric redshift distribution (peak of the probability function for each source) for the selected excess sources for 
COSMOS and UDS. Both distributions peak at z ~ 0.85, as expected, thus implying that the great majority of the selected sources are Ho; emitters. The right 
panel shows the best photometric redshifts plotted against the spectroscopic redshifts for the emitters with accurate redshifts (reliability > 95%) measured in 
the 2;-COSMOS DR2 dataset. Dash-dotted lines indicate the redshifts at which emission lines are matched with the NB j filter and these are the same in both 
figures (the width of the lines scales with the fraction of confirmed emitters; see Table|2]for more details). The great majority of potential emitters is confirmed 
to be real and only a small fraction (~ 6%) was not confirmed. 



Table 2. Spectroscopic redshifts matches (from 2;-COSMOS DR2) for the narrow-band excess objects. The first column of the table indicates the emission line 
producing the narrow-band excess, based on the accurate spectroscopic redshifts which place such lines within the narrow-band filter. The remaining columns 
present the mean redshift, the number of emitters and the number of those selected as Ho; using the photometric redshift selection (see §3.4). 
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between the selected Ha emitters and the remaining emitters, re- 
vealing both a population of lower and higher redshift. The right 
panel of Figure |4] also reveals that the Ha selected emitters present 
colours which clearly distinguish them from the low-z and higher-z 
samples of emitters - they also occupy a region in the colour-colour 
diagram where one expects to find star-forming galaxies at z=0.84. 
Finally, as examined below, spectroscopic redshifts for ~ 25% of 
the sample in COSMOS show that only a negligible fraction of the 
Ha emitters are being missed by the photometric redshift selection. 



3.5 Spectroscopic redshifts and selection robustness 

Both COSMOS and UDS have large spectroscopic surveys under- 
way, but only a limited fraction of those spectroscopic redshifts is 
currently available. Nevertheless, for the ^-COSMOS survey ( Lilly 
let al.pOO?] ) Data Release 2 (DR2), the match with the sources 
detected on the NBj frames yields 4600 sources (using a 1.0^^ 
match criteria), of which 138 (each matching only one z-COSMOS 
source) have been selected as excess sources and hence potential 



emitters. Twelve of these have unreliable or non-existent redshift 
determinations. Of the remaining 126, the vast majority (119) have 
redshifts which place an emission line within the NB j filter (see Ta- 
ble [2] for the detailed list of emitters), confirming the narrow-band 
excess as an emission line. For the other 7 sources, it was not possi- 
ble to identify any emission line falling into the NB j filter. Spectra 
for these sources were analysed in order to look for potential errors 
in the redshift determination. In fact, while redshifts for 5 of them 
are very robust, one source, identified as 2; ^ 0.4, seems better fit- 
ted by being at z = 0.85 (with [Oil] 3727 being detected, while 
the fit at lower redshift assumes an emission line where fringing 
starts to become an issue in the spectrum). The remaining source 
contains one emission line (identified as [Oil] 3727) but with a low 
signal to noise (S/N < 2.0). These 7 sources all have a colour ex- 
cess significance S = 2.5-3.0, and at least the 5 robust cases repre- 
sent the galaxies expected to be randomly scattered into the candi- 
date emitter list at these low E values. All of the candidate emitters 
with E > 3.0 were confirmed to be real. It is also worth notic- 
ing that ^-COSMOS DR2 is highly biased towards z < 1 sources, 
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Figure 4. The left panel shows the photometric redshift distribution (peak of the probability function for each source) for the selected Ho; emitters (see Section 
\3A] and the remaining emitters, with the spectroscopically confirmed lines plotted as in the left panel of Figure [3] The right panel presents a BRIK colour- 
colour diagram, clearly distinguishing the high-redshift emitters from the low-redshift emitters, and also the Ho; emitters from these two populations. Overlaid 
on the plot are evolutionary tracks for different galaxy models, drawn from Poggiantij|1997|; small symbols on the tracks mark steps of 0.5 in redshift from 
to 2.5, whilst the larger symbol indicates z=0.84. The observed colours for the H-alpha emitters are consistent with them being mostly star-forming galaxies 
at z=0.84, whilst the colours of the other emitters also lie in the expected regions. 



which means that its completeness for potential interlopers (which 
are more likely to be in this redshift range) is very high, while at the 
same time it misses most of the emitters at higher redshift. There- 
fore, the contamination rate within the sample of emitters (i.e., the 
fraction of non-emitters) is likely to be lower than ~ 6%, even 
down to S = 2.5. 

Tests regarding the selection criteria used were also done by 
selecting samples with different colour significances and different 
equivalent width cuts, and comparing the matches for each selec- 
tion criteria with z-COSMOS DR2. As expected, the number of po- 
tential interlopers declines rapidly with an increasing colour signif- 
icance threshold, but so does the sample size. The sample with the 
largest number of real emitters whilst returning a low ratio between 
potential interlopers and real emitters is obtained using E 2.5 
and an equivalent width cut of ~ 50A, therefore indicating that the 
selection criteria used are producing reliable results. 

Table [2] compares the spectroscopic selection of line emitters 
with the selection from the photometric redshift catalogue. It can be 
seen that 88 of the 93 Ha emitters are correctly identified 95% 
completeness). On the other hand, 2 non-Ha emitters (with fluxes 
which do not deviate significantly from the median) were found in 
the photometric redshift selected sample (implying 3. ^ 98% reli- 
ability), with one of those being a [Sll] 6731 emitter. Because this 
line is very close to Ha it is very difficult to completely distinguish 
these emitters with photometric redshifts alone. However, as seen 
by the sample which has been presented, the total contamination by 
these will be very small (< 2%). The spectroscopic data were used 
to improve the Ha catalogue, removing the 2 non-Ha emitters and 
including those 5 which had not been selected by the photomet- 
ric redshifts. These few sources present mean and median fluxes 
slightly higher than the rest of the sample, suggesting that some 
bright Ha emitters could be lost by a photometric redshift selec- 
tion, due to the emission line contribution to the broad-band fluxes. 
However, as the brighter Ha emitters are likely to contain strong 
emission lines in the visible ([Oil] 3727, H/3 and [Olll] 5007), these 
targets should have a very high completeness in z-COSMOS, so 



this can be interpreted as an upper limit. As these emitters are in- 
troduced into the sample, no further specific correction was applied. 

While 93 Ha emitters are confirmed using z-COSMOS, there 
are 72 z-COSMOS galaxies with a spectroscopic redshift which 
should place the Ha line in the NBj filter. The failure to se- 
lect these as narrow-band excess sources means that they have a 
weak or absent Ha emission line. These present a mean colour 
(J— NBj) = 0.1 corresponding to a mean measured Ha flux 
lower than the survey limit (3 x 10~^^ erg s~^ cm~^ and a mean 
EW of 20A). The brightest source has J = 19 mag. These are 
likely to be very faint emitters for which a completeness correc- 
tion is applied later (see §4.1. 2| ). The [Oil] 3727 line fluxes from 
z-COSMOS galaxies confirms this: where they are detected they 
yield similar line fluxes to the ones that would be estimated for 
Ha. On the other hand, some of these sources have a negative 
(J— NBj) colour. By applying similar selection criteria used for 
emitters, it seems that at least 2 Ha absorbers would be selected. 



3.6 Narrowband-iT and H matches 

By design, the custom-built narrow-band filters make it possible 
to look for line detections in multiple bands to refine redshift esti- 
mates. Currently, data are available in H2S1 for COSMOS (Geach 
et al. 2008) and NBh and H2S1 for UDS (Sobral et al., Geach et 
al., in prep.) 

For COSMOS, it was possible to match 3 sources between 
excess sources catalogs (J and K), confirming a z = 2.23 redshift 
for those emitters. This shows that the line being detected in NB j 
is [On] 3727. 

For UDS, 18 sources were matched between the NBh and 
NBj. These have photometric redshifts consistent with either be- 
ing Ha emitters at 2; ^ 1.5 or Ha emitters at 2; = 2.23, thus 
meaning that the line detected in the J band is H^^ or [Oil] 3727 
(respectively) for these emitters. At least 2 of these are also selected 
as emitters in the NBk filter, which confirms z — 2.23 for them. 
In addition, there are another 6 matches between the K and the J 
narrow-band filters, indicating that those sources are at 2; = 2.23 
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as well. Four of these 24 matched sources were found in the pho- 
tometric redshift-derived Ha catalog and were thus removed as the 
evidence points towards these being at higher redshift. 



4 RESULTS 

A total of 743 candidate z — 0.84 Ha emitters are detected 
over 1.30 deg^ down to an average observed line flux limit of 
8 X 10~^^ergs~^ cm~^. This sample will now be used to eval- 
uate the Ha luminosity function and estimate the star formation 
rate density at z = 0.84. The morphological mix of these Ha emit- 
ters in the COSMOS field will then be investigated, together with 
their evolution as a function of luminosity and redshift, and their 
contribution to the total luminosity function and star formation rate 
density. 



4.1 Ha luminosity function at z = 0.84 

In order to calculate the luminosity function of Ha emitters, line 
fluxes are converted to luminosities by applying: 



100 



(6) 



where Dl is the luminosity distance, 5367 Mpc at z = 0.84. 

The estimate of the source density in a luminosity bin of width 
A(logL) centered on logLc is given by the sum of the inverse 
volumes of all the sources in that bin. Therefore, the value of the 
source density in that bin is 



0(log(ic)) 



A(logL) 



E 



:l< 



A(logI.) 



A{Vmte 



(J) 



Here, i refers to sources and c to the center of each bin. The vol- 
ume probed is calculated taking into account the survey area and 
the narrow-band filter width (initially assumed to be a top-hat func- 
tion across 1. 2037-1. 2185^Lm; see §4.1.4| for a refined approach), 
resulting in a co-moving volume of 7.4 x 10^ Mpc^ (UDS) and 
1.04 X lO^Mpc^ (COSMOS). As detailed before, this takes into 
account the removed area due to the presence of bright stars and 
consequent artifacts, along with noisy areas. 

The luminosity functions presented here are fitted with a 
Schechter function defined by the three parameters: a, 0* and L* : 

(i){L)dL = (i)\L/L''Yeyi^{-L/L'')d{L/L''). (8) 



4.1.1 [ Nil] flux contamination correction 

When computing line fluxes and equivalent widths for the Ha line, 
one must note that the adjacent [Nil] lines at 6548 and 65 83 A will 
also contribute to both quantities, increasing them both. In order 
to account for this, the relation between the flux ratio F[nii] /Fh^ 
and the total measured equivalent width EW(Ha-i-[Nll]) derived by 
Villar et al. (2008) is used. This shows that the fractional [Nil] flux 
contribution decreases with increasing EW. The relation between 
y = log(F[Nii]/FHa) and x = log(EW(Ha + [Nil])) can be 
approximatedby ^ = -5.78 + 7.63x- 3.37x2 + 0.42x^ forEW> 
50A. This results in a correction lower than the conservative 33% 
used by some authors (e.g. Geach et al. 2008), with the median 
being 25% for this sample with EW > 50 A. 
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Figure 5. The total completeness function used for COSMOS and UDS 
resulting from several simulations to address the effects of selection and 
detection. This confirms a ~ 30% completeness at a flux limit of ~ 
10-1^ erg s-^cm-^. 



4.1.2 Completeness correction 

Fainter sources and those with weak emission lines might be missed 
and thus not included in the sample; this will result in the underesti- 
mation of the number of emitters, especially at lower luminosities. 
Furthermore, the completeness rate is highly affected by the selec- 
tion criteria (E > 2.5, EW > 50A, cr > 3). To address this prob- 
lem, a series of simulations were conducted. First, the recovery rate 
has been studied as a function of magnitude for objects in the same 
frame, using a Monte Carlo method. For this, 10 different galax- 
ies (both real - taken from the narrow-band images in UDS and 
COSMOS - and simulated) were used. These were introduced (20 
for each run, in any given image) with different input magnitudes 
into the science frames and then detected using the same extrac- 
tion parameters as the main catalogues. The recovery rate and the 
recovered magnitude were then studied. The latter follows the in- 
put magnitude reasonably well down to the 3-a limit. The recovery 
rate, while varying slightly with the type of galaxy used (the arti- 
ficial and point-like objects, for example, showed a much higher 
recovery rate at fainter levels than real galaxies), falls off sharply 
fainter than NB j ~ 21. Similar simulations were also done for the 
J band images. 

This information alone cannot be used to correct the luminos- 
ity function, as faint narrow-band sources don't necessarily have 
faint Ha lines and vice-versa. To address this, a second series 
of simulations were performed to investigate the selection com- 
pleteness of faint emission lines within detected galaxies. For this, 
narrow-band detections not classified as emitters were considered 
not to have any emission line. For each of those sources 21000 
for COSMOS and - 15000 for UDS) a line flux (in steps of 
5 X 10~^^ erg s~^ cm~^ up to 5 x 10~^^ erg s~^ cm~^) was added 
to both NB J and J-band magnitudes, and the selection criteria were 
applied to the revised magnitudes for each line flux to study the re- 
covery rate. This outputs a colour selection completeness for each 
line flux. Finally, a combined completeness correction was done 
combining both simulations described. The combined complete- 
ness function, which is used to correct the luminosity function on a 
bin to bin basis, is shown in Figure[5] Slightly higher completeness 
rates at the same flux level have been claimed in | Villar et al.| ( |2008"] ), 
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for shallower data, but it should be pointed out that these authors 
use a combination of selection criteria which is significantly dif- 
ferent from those used in this study. On the other hand, they only 
consider the colour- selection completeness, neglecting the detec- 
tion completeness, and on the other hand, they use sources detected 
down to just 0.8a in NBj (instead of the 3-a limit used in this 
work), and with a lower EW 15 A) cut; therefore, the results 
presented here are likely to be more reliable than those by |Villar| 
let al., ( ,2008| ). 

4.1.3 Extinction Correction 

It is well-known that the Ha emission line is not immune to dust 
extinction, although it is considerably less affected than Lya or the 
UV continuum. Measuring the extinction for each source can in 
principle be done by several methods, ranging from spectroscopic 
analysis of B aimer decrements to a comparison between Ha and 
far-infrared determined SFRs. For now, however, a conservative 
^Ha = 1 mag is used, which is the same correction adopted in 
most of the similar studies done before (e.g. jPujita et al.|2()03{[I^ 
|et al.|[2007] |Geach et al.||2008| ). This same correction is also ap- 
plied to all other datasets being compared here, when possible, in 
order to give consistency between all studies and search for evolu- 
tion. §4.2.1| shows that this correction should not be too far from 
the actual mean extinction for the entire sample. 



4.1.4 Filter profile corrections 

The narrow-band filter transmission function is not a perfect top- 
hat (as assumed earlier), and thus the real volume probed varies 
as a function of intrinsic luminosity: luminous Ha emitters will be 
detected over a larger volume than the fainter ones because they 
can be detected in the wings of the filters (although they will be 
detected as fainter sources in these cases). Low luminosity sources 
will only be detected in the central regions of the filter and thus the 
effective volume will be smaller. 

In order to correct for this when deriving the Ha luminosity 
function, a further set of simulations was run. Firstly, the luminos- 
ity function was computed with the corrections described above 
and the best fit was assigned. This was then assumed to be the true 
luminosity function, allowing the generation of a set of ~ 10^ Ha 
emitters with a flux distribution given by the measured luminosity 
function, but spread evenly over the redshift range z = 0.81-0.87 
(assuming no cosmic structure variation or evolution of the lumi- 
nosity function over this narrow redshift range). The top-hat filter 
model was then confirmed to recover the input luminosity func- 
tion perfectly. Next, the true filter profiles were used to study the 
recovered luminosity function. These simulations showed that the 
number of brighter sources is underestimated relative to the fainter 
sources. A mean correction factor between the input luminosity 
function and the one recovered (as a function of luminosity) was 
then used to correct each bin. The simulation was run again with 
the new luminosity function, confirming that the recovered lumi- 
nosity function is very similar to the input luminosity function. 

The filter profiles were also checked against the spectroscopic 
redshift distribution from z-COSMOS DR2. By distributing the ar- 
tificial Ha emitters to reproduce the luminosity function, one can 
predict the redshift distribution of observed sources for each filter. 
Those matched well the redshift distribution of the 2;-C0SM0S 
emission line sample, although this sample only contains ^ 100 
emitters at the moment. 



4.1.5 Fully corrected luminosity function 

The final z — 0.84 Ha luminosity function is presented in Figure|6] 
and in TablejS] The raw luminosity function, without correcting for 
incompleteness and filter profile biases (correcting only for [Nil] 
contamination and Aho. — 1 mag of extinction), is also shown. 
The errors are Poissonian in each bin, combined with an uncertainty 
on the correction factor assumed to be 10% of the applied correc- 
tion. Luminosity functions are also computed separately for the two 
observed fields (Figure|6]) with these being compared with the com- 
bined luminosity function. Both luminosity functions are generally 
consistent within the errors, although there seems to be a slightly 
higher density at z = 0.84 in COSMOS. The on-sky distribution of 
the Ha emitters in the COSMOS and UDS fields is far from homo- 
geneous: they are highly clustered and several high-density regions 
can be found (a detailed clustering analysis is underway and will 
be published in Sobral et al. in preparation). Assuming a standard 
angular correlation function for the Ha emitters, parameterised as 
w{0) = A0~^-^ (with in arcsec), then for a value of A ^ 1 
(typical for galaxies like these), and following [Peebles | ( | 1 975 1 ) , the 
cosmic variance in the number of emitters expected in a 1.3 deg^ 
sky area due to clustering is estimated to be roughly double the 
pure Poissonian errors. This leads to (lcr)uncertainties of ^ 12% 
on the total number of detected Ha emitters; the corresponding un- 
certainties for the individual fields are ^15% for UDS and ^\3% 
for COSMOS. These are fully consistent with the actual differences 
in number densities found. Even for an extremely clustered popu- 
lation, with A 10 — 15, the source count variations would not be 
larger than ^ 20% on fields of this size. Given this, and the consis- 
tency between the COSMOS and UDS fields (and also the results 
of Villar et al 2008), it is safe to say that the results derived are 
largely robust against cosmic variance. 

Also plotted in Figure [6] are other published Ha luminosity 
functions at different redshift; these demonstrate a significant evo- 
lution with redshift, linked to an increase of the number density of 
both faint and bright emitters at least up to z 1. To quantify this 
evolution a Schechter function is fitted to the combined z = 0.84 
luminosity function, which yields: 0* = 10~^-^^^°-^° Mpc~^, 
a = -1.65 ±0.15 and L* = 10^^-^^^° °^ erg s-\ These param- 
eters, and a comparison with the best fits for COSMOS and UDS 
separately, can be found in Table 5. The best fit parameters indicate 
a steeper faint end slope at z = 0.84 than the canonical a = —1.35 
usually assumed; similarly steep faint-end slopes for the luminos- 
ity function of star-forming galaxies have recently been found from 
UV studies (e.g. |Bouwens et al. 2008). These authors have argued 
for an increase in the faint-end slope of star-forming galaxies with 
redshift. However, evidence for a steep faint-end slope in the Ha 
luminosity function (a ^ —1.6) has been found even at z = 0.08 
by Ly et al. (2007) as well as in other Ha studies at higher red- 
shifts, z - 0.2-1.3 (e.g . Hopkins et al.||2000l [Fi^ta et al. 2003^ 
[Dale et al.|2008| ). Figure [7] shows the best-fit a for Ha luminosity 
functions derived at ^ = 0.02-1.3. This reveals no clear evidence 
of evolution of the faint-end slope with cosmic time; indeed, rather 
than having a random distribution, the values of a seem to clus- 
ter into two groups with a ^ —1.6 or a ^ —1.3. The origin of 
these discrepant results is unclear, but cosmic variance may play an 
important role. 

Figure [7] also shows the variation of L* and 0* of the Ha lu- 
minosity function out to z ^ 2. To derive these values, only LPs 
of Ha surveys which were fitted with the "canonical" a=— 1.35 
or that had their data points published (those were then refitted by 
fixing a=— 1.35) were used, in order to reduce degeneracies and 
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Table 3. The Ho; luminosity function. Superscript numbers in the column titles indicate the corrections made in each stated value: 1 - extinction correction 
{^Ha. = 1 niag) and [Nil] correction as a function of total measured EW; 2 - completeness correction; 3 - filter profiles biases correction. 
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Figure 6. The left panel shows the Ho; luminosity function (corrected for [Nil] contamination, completeness, extinction and filter profile biases) with the best 
fit Schechter function given by 0* = io-i-92±o.io Mpc-^, q, = -1.65 ± 0.15 and L* = io42.26±o.05 ergs-^ Small black dots present the derived 
luminosity function with corrections for [Nil] and extinction {Ah a. = 1 mag) only. Other luminosity functions (corrected for Ana = 1. except for z = 0.08 
and z = 0.24 where the extinction corrections from the authors were used) from narrow-band surveys of Ho; emitters at different redshifts are presented 
for comparison, showing a clear evolution with cosmic time up to at least z ^ 1. The derived Ho; luminosity function for each separate field (COSMOS 
and UKIDSS UDS) is presented in the right-hand panel, with a comparison with data points from Villar et al. (2008) - assuming the same dust extinction as 
applied in the present study. There seems to be a slightly higher density of sources in COSMOS, but both luminosity functions agree well with the combined 
fit (dashed line). 
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Figure 7. The evolution of the best-fit Schechter function parameters of the Ho; luminosity function for Ho; surveys. The left panel presents the best-fit values 
for the faint end slope, a, which appear to fall into two distinct groups with a ~ —1.6 and a ~ —1.3, with little evidence for any redshift evolution. The 
middle and right-hand panels show the evolution of L* and as calculated by fixing a = —1.35. L* appears to evolve strongly with redshift out to z > 2, 
whereas peaks at z ~ 1 and then declines. Fixing the faint-end slope at o; = —1.65 yields a similar result. 
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allow a direct comparison between results. It should be noted that 
some degeneracy will still remain between L* and 0* values. Fig- 
ure [7] shows a strong evolution in L*, increasing by at least an or- 
der of magnitude from the local Universe (Gallego et al. 1995) to 
z = 2.23 (Geach et al. 2008). The evolution of 0* is somewhat 
different: while it appears to increase from z = up to 2; = 0.84 
(this work and Villar et al. 2008) by one order of magnitude, it 
would then need to fall at 2; > 1 to be consistent with the higher 
redshift data of Yan et al. (1999) and Geach et al. (2008). Note that 
the increased value of 0* at z = 0.84 does not arise just due to the 
degeneracy between 0* and L*: a value of 0* = 10~^"^ Mpc~^ 
for the current data- set can be strongly rejected, having a proba- 
bility < 10~^. The evolution of L* and 0* is therefore revealing 
important details of the evolution of the Ha luminosity function: 
from 2; ^ to 2; 1 it seems to be driven by an increase in the 
number density of both bright and faint emitters, with an increasing 
population of bright emitters at higher redshift then being responsi- 
ble for a stronger evolution in L*. These results are also consistent 
with studies done using 24 /xm data (e.g. |Caputi et al.|2007| ). The 
evolutionary trends do not change if a is fixed at a higher value 
(a = —1.65 for example). 

4.2 The star formation rate density at z = 0.84 

4.2. 1 AGN contamination 

The Ha luminosity function previously derived used all the Ha 
emitters from the survey, and while most of such sources are likely 
to be star-forming galaxies, some of these can also be AGN. Spec- 
tra from 2;-C0SM0S DR2 were used to explored this. A visual in- 
spection of the 93 available spectra was done to confirm additional 
emission lines ([Oil] 3727, [Olll] 5007 and H^^) and the assigned 
redshift. The line fluxes were then measured using an IDL script. 
Although the comparison of those lines with the Ha line fluxes is 
influenced by many factors (e.g. Ha/[Nll] ratio, exact location of 
the Ha line within the filter profile, fraction of emission line light 
falling into 2;-C0SM0S slit) it is noteworthy that both the mean 
ratio of Ha/[Oll] 3727 = 2.27 and the ratio of Ha/H/3 = 4.16 are 
consistent with an Ha extinction of ^ 1 mag or slightly higher. 

In order to estimate the AGN contamination, the 
[On] 3727/H;5 and [Olll] 5007/H;^ line ratios were used; these 
have been widely used to separate AGN from star-forming galaxies 
(e.g. |Rola et al.|1997] ). Only spectra with all lines being detected 
at S/N > 3.0 were used, which results in a sample of 28 galaxies, 
mainly due to the low S/N at longer wavelengths where [Olll] 5007 
and H/3 are found. Figure [8] shows data-points for the line ratios, 
while the curves represent maximum line ratios for a star-forming 
galaxy (from OB stars with effective temperatures of 60000 K and 
50000 K). From the sample of 28 Ha emitters, 23 seem to be clear 
star-forming galaxies, while 3 are likely to be AGN contaminants. 
A ^ 15% AGN contamination is thus estimated, consistent with 
that found in other Ha studies. The AGN are found to have Ha 
fluxes typical of the rest of the sample. 

4.2.2 Star formation rate density 

The observed Ha luminosity function can be used to estimate the 
average star formation rate density, psfr.^^z = 0.84. To do this, 
the standard calibration of Kennicutt (1998) is used to convert the 
extinction-corrected Ha luminosity to a star formation rate: 

SFR(M0year"^) = 7.9 x 10"^^ Lhc (ergs"^). (9) 
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Figure 8. Line ratios from the z-COSMOS spectra of the z = 0.845 Ha 
sample. These show that the great majority of the sample is composed of 
star-forming galaxies (82%), as expected, with 11% showing evidence for 
being AGN contaminants and 7% being unclassified. The red curves repre- 
sent the maximum line ratios for a star-forming galaxy (from OB stars with 
effective temperatures of 60000 K (solid line) and 50000 K (dashed line)). 

This assumes continuous star formation. Case B recombination at 
Te = 10^ K and a Salpeter initial mass function ranging from 
0. 1-100 M©. All measurements of psFR include a correction of 
15% for AGN contamination and an extinction correction Ah^ = 
1 mag, except where the authors only presented their own extinc- 
tion corrected luminosity function. 

In §4.1.5| a significant evolution in the observed Ha lumi- 
nosity function was observed. The left panel of Figure |9] shows 
how this translates into an evolution in ps fr as a function of red- 
shift, for luminosity functions which have been integrated down to 
hua > lO^^^ergs"^ (the limit of this survey). The measurement 
at 2; = 0.84 presented in this study (0.15 ± 0.02 M© yr~^ Mpc"^) 
demonstrates a strong rise in psFR, when compared to the local 
Universe { [Gallego et al.||1995| |Perez-Gonzalez et al.||2003[ |Ly et 
"ah 2007) and low redshift measurements (e.g. Tresse & Maddox 
.1998, Sullivan et al. 12001 1 [Dale et al.|2008||Morioka et al. 2008^ 
|Westra& Jonesr 2008; Shi oya et al.|2008l|Sumiyoshi et al.|200"9l ), 
as suggested by other smaller surveys done at similar redshifts (e.g. 
Tresse^ etaLl2002l [Villar et al.|2 008). This rise seems to be shghtly 
steeper than Psfr {'^ + z)^. When compared to higher redshift 
(e.g. Geach et al. 2008), the observations also support a flattening 
in psFR around 2; ~ 1, up to at least z = 2.23. A rise and subse- 
quent flattening of the star formation rate density out to 2; '-^ 2 has 
therefore been accurately measured using a single star formation 
tracer. Cosmic evolution of dust reddening corrections may alter 
the results slightly but would have to be very strong to change the 
overall conclusions. 

Figure [9] also presents the same evolution, but now integrating 
the entire luminosity function. Caution should be used in interpret- 
ing this figure as it involves extrapolating all the luminosity func- 
tions and it is critically dependent on the assumed faint-end slope. 
For this study, for example, psfr = 0.37 ± 0.18 M© yr~^ Mpc~^ 
using the measured value of a = —1.65, but if one adopts a = 
-1.35, PSFR is reduced to O.28±O.O8M0 yr~^ Mpc~^. Measure- 
ments obtained from a set of other star-formation indicators com- 
piled by Hopkins (2004) are also shown for comparison (corrected 
by a common extinction factor consistent with the Ha extinction 
correction applied here). This confirms the same rise seen using 
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Figure 9. The evolution of the star formation rate density as a function of redshift based on Ho; - down to the limit of this survey (left panel) and by 
integrating the entire luminosity function (right panel). The right panel also includes data from all other star-formation tracers, as detailed in the Hopkins 
(2004) compilation. Both panels show a clear rise in psFR up to at least z ~ 0.9-1.0, slightly steeper than the canonical (1 + z)^ rise, followed by a 
flattening out to at least z ~ 2.2. Ha data points refer to (in order of increasing mean redshift) Gallego et al. (1995; z = 0.02), Perez-Gonzalez et al. (2003; 
z = 0.02), Ly et al. (2007; z = 0.08), Sullivan et al. (2000; z = 0.15), Dale et al. (2008; z = 0.16), Tresse & Maddox (1998; z = 0.20), Shioya et al. (2008; 
z = 0.24), Fujita et al. (2003; z = 0.24), Dale et al. (2008; z = 0.24), Morioka et al. (2008; z = 0.24), Ly et al. (2007; z = 0.24), Ly et al. (2007; z = 0.4), 
Tresse et al. (2002; z = 0.7), Villar et al. (2008; z = 0.84, 1 mag extinction - higher psFR - and author's extinction correction), this study (z = 0.845), 
Hopkins et al. (2000; z = 1.25), Yan et al. (1999; z = 1.3), Geach et al. (2008; z = 2.23). 



only Ha and the possible flattening at z ^ 1. However, it appears 
the Ha measurements at z > 0.7 are finding slightly higher values 
of ps FR when compared to the average of all other measurements, 
while at z < 0.3 they measure slightly lower values. This may re- 
flect a number of systematic errors or factors such as evolution in 
the typical reddening in the star forming population. 



4.3 The morphology of Ha emitters 

4. 3. 1 Morphology with HST imaging 

The COSMOS field has sensitive Hubble Space Telescope (HST) 
ACS F814W coverage which provides detailed morphological in- 
formation on the Ka sample. These data are used to study the mor- 
phologies of the selected star-forming galaxies. This morpholog- 
ical information can be greatly enriched with colour information, 
and thus the ACS images are combined with deep Subaru data in 
Brz bands to produce a pseudo-true colour image. The Subaru im- 
ages are first registered and transformed to match the HST images 
using IRAF and PYTHON scripts. The HST image is then used to 
define the luminance of the true colour image. In this way colour 
and morphological information were derived on ~ 500 Ha emit- 
ters with sufficient detail that a visual morphology analysis could 
be undertaken. Figure [TO|presents some examples of these. 

All of the galaxies within the sample were visually classified. 
In order to compare those results with morphologies obtained in 
an automated way (by ZEST, kindly supplied by Claudia Scarlata, 
[Scarlata et al.| ( |2007] )), 4 morphological classes were used: 1) Early- 
types, 2) Disks/spirals, 3) Irregulars, and 4) Unclassified. Table [4] 
shows the comparison between ZEST and the visual classification, 
while Figure [To| shows examples of each morphological type. In 
general, a very good agreement between the visual classifications 
and ZEST was obtained. The main difference is that visual clas- 
sification was able to reduce the number of unclassified galaxies. 
Also, visually, it was possible to improve the distinction between 
disks and irregulars by using the colour information - distinguish- 



ing between multiple bulges and bright spiral arms - which resulted 
in a sightly lower fraction of irregulars in the visual classification. 

Due to a considerable number of sources showing evidence 
of merging activity, the sample was also classified independently 
into merger classes. A source is classified as a merger when: 1) it 
presents a clearly disturbed morphology or disturbed disk which 
is inconsistent with being that of a normal disk galaxy; or 2) 
the source presents more than 1 bright point-like source and the 
colour information is inconsistent with one of those being a spi- 
ral arm; or 3) there are two or more galaxies which are very close 
(< 15kpc). By applying these criteria, a visual classification was 
done, where each source was classified at least twice. In the end, 
3 merger classes where used to classify all sources: the consistent 
non-mergers (non-mergers); those galaxies classified at least once 
as probable mergers (potential mergers); and galaxies always clas- 
sified as mergers (mergers). The distribution of mergers within the 
previous morphological classifications and the total numbers can 
be found in Table|4]and examples of classified galaxies can also be 
found in Figure fToT 

From a total of 477 Ha emitters, 381 are disks (80%), with 68 
being classified as irregulars (14%), 14 are early-types (3%) and a 
total of 13 (3%) are unclassified as these are too faint. Furthermore, 
24% of the sample seems to be populated by clear mergers, with 
the total merger fraction being estimated as 28 =b 4%; almost all 
irregulars fall into this class. 

4.3.2 Morphology-Ha luminosity relations 

The left panel of Figure [TT] presents the fraction of galaxies classi- 
fied into each morphology as a function of Ha luminosity. There is 
a clear evolution of morphological type with increasing Ha lumi- 
nosity. While at low star-formation rates, disks dominate the sample 
completely (> 85%), at higher Ha luminosities (L > L*) irregu- 
lars become more significant, reaching 100% in the highest lumi- 
nosity bin. With the large sample presented in this work it is pos- 
sible to derive independent luminosity functions for each of those 
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Table 4. Visual morphology classifications compared to the automated classifications obtained by ZEST (Scarlata et al. 2007) for the Ho; emitters at z=0.84. 
The two agree very well. Differences arise mostly from using colour information for the visual classifications, while ZEST only uses the HST F814W imaging. 
The final rows of the table give the distribution of visual morphologies for different merger classes. 
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Class 


Early 


Disks 
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Unclassified 
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Non-mergers 
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Early 


7 


3 





4 


14 (3%) 


12 





2 


Disks 


2 


301 


15 


63 


381 (80%) 


297 


33 


51 


Irregulars 





10 


44 


14 


68 (14%) 


1 


7 


60 


Unclassified 





2 





12 


14 (3%) 











Total 


9 (2%) 


316(67%) 


59 (12%) 


93 (19%) 


477 (100%) 


310(67%) 


40 (9%) 


113(24%) 



morphological classes. Those can be seen in Figure [12] with the 
best fit Schechter function parameters tabulated in Table 5. They 
illustrate well the different contributions to the total Ha luminosity 
function by disks and irregulars. Irregulars present a remarkably flat 
luminosity function, only falling at the highest luminosities, while 
disks demonstrate a steeper faint end and a much lower space den- 
sity at the bright end. Overall disks are the dominant contributors 
to the total psfr, with a large number of galaxies producing stars 
at rates < lOM0yr~\ 

As mentioned before, almost all of the irregulars show clear 
evidence of merging activity. In fact, separating non-mergers from 
likely mergers (mergers, together with potential mergers weighted 
by 0.5), a very similar behaviour is found (Figure pT|. While non- 
mergers dominate at faint luminosities, the contributions from the 
two populations cross over at L ~ L*, with mergers dominating 
the bright-end of the luminosity function. The luminosity functions 
for non-mergers and potential mergers are presented in Figure [12] 
This reveals that non-mergers present the typical disk luminosity 
function found before, but mergers have a steeper luminosity func- 



tion than the irregulars. This means that while irregulars on their 
own are only important in the bright end of the LF, mergers seem 
to play a dominant role at the bright end (accounting for^^ 60% of 
the psFR there) together with a non-negligible contribution even 
at the faintest luminosities. Overall, mergers account for ^ 20% of 
the total psFR z — 0.84. 

The relations which have been presented were found to be 
very robust, with the same results being obtained regardless of the 
use of the visual or ZEST classifications, and independently of the 
[Nil] correction - a constant correction also produces the same de- 
pendences. 

433 Redshift evolution of the morphology relations 

One can compare the morphological mix for the Ha sample at 
z — 0.84 with that found at lower redshifts. In the local Universe, 
the morphologies of the Ha emitters seem to be somewhat differ- 
ent, with visual morphologies from |Vitores et aL| ( |1996| ) indicating 
that disks are 88% of the sample, while early-type galaxies increase 
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Figure 11. Morphological class fractions (left panel) and merger fractions (right panel) as a function of Ho; luminosity at z = 0.84. These show clear Ho; 
luminosity dependences, with the fraction of irregulars and mergers rising with increasing luminosity. The vertical line indicates the value of L* derived from 
the best total Schechter function fit with a = —1.65. 




Figure 12. The computed luminosity functions and the best fits for different visual morphologies. The dash-dotted line refers to the z = 0.08 LF (Ly et al. 
2007), with the solid Hne being that at z = 0.24 (Shioya et al. 2008), and the dashed line is the best fit for the combined UDS+COSMOS luminosity function 
2it z = 0.84. The brightest bin for both irregulars and potential mergers seems much higher than in the combined luminosity function due to the COSMOS 
field having a higher space density of bright sources (c.f. Figure 5). 



their importance (9%) and irregulars drop to only 3%. At z = 0.24 
the narrow-band Ha emitters of Shioya et al. (2008) in the COS- 
MOS field can be more directly compared to the 2; = 0.84 sam- 
ple, using ZEST morphologies for both (although noting that the 
rest-frame wavelength of the images is different between the two 
studies - rest-frame B band at 2; = 0.84 and rest-frame R band 
at 2; = 0.24 - which could introduce a small bias). In the ZEST 
classifications at z = 0.24, disks are still dominant (89%), ^-^10% 
of the galaxies have an irregular morphology, and elliptical galax- 
ies account for only 1 % of the sample. These results seem to show 
that while disks are always the dominant population (> 80%) up to 
z ^ 1, there is a significant increase in the irregular fraction from 
3% at z = to 10% at z = 0.24 and 15% at z = 0.84. 

The ZEST classifications for the Shioya et al. (2008) sample 
can also be compared with those for the 2; = 0.84 sample to see 
if there is any evolution in the morphology-Ha luminosity rela- 
tions. Thus Figure \T3\ presents the comparison between z = 0.24 



and z = 0.84 with the luminosities scaled by L*. A very simi- 
lar variation in morphological mix with Ha luminosity is found at 
z — 0.24, with the irregular fraction increasing from ^ 10% to 
^35% with increasing luminosity and with the disk fraction hav- 
ing the opposite behaviour, decreasing from ^ 90% to ^ 60%. The 
relatively good agreement between the data at these two different 
redshifts and the fact that the switch-over luminosity appears to be 
at the same L/L* at both redshifts seems to point towards a rather 
simple L* evolution driving the morphology-Ha luminosity rela- 
tion. However, the Shioya et al. (2008) sample is not able to probe 
the brightest sources and thus there is no low redshift data to di- 
rectly compare the complete dominance of irregulars at L > 3L*. 

The disk fraction at z = 0.24 appears to be slightly lower 
than the disk fraction at 2; = 0.84 (for the bins which can be di- 
rectly compared). Curiously, this is due to an apparent rise in the 
fraction of early-types. However, this might be simply reflecting 
classification errors between the two epochs. At least some of the 
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Table 5. The Ho; Luminosity function and psFR for different subsets of the data. Both \ogpL >4i.5 and psFR >4i.5 result from integrating the Ha 
luminosity function down to 10^^ -^ ergs~^. psFR total results from a complete integration of the luminosity function. A 15% AGN correction is made to 
PSFR measurements. 
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Figure 13. Morphology class fractions as a function of Ho; luminosity for 
z = 0.24 (Shioya et al. 2008) and this study. These results find the same 
morphological dependence at both epochs consistent with an evolution in 
L* (log(L) ~ 0.4). Morphologies were obtained in an automated way by 
ZEST (Scarlata et al. 2007) for both samples. 



very faint galaxies which have been unclassified at z = 0.84 (even 
after visual classification) are likely to be early-types which are just 
too faint to be seen at 2; = 0.84, but sufficiently bright dX z — 0.24 
to be classified. 



5 CONCLUSIONS 

Deep near-infrared narrow-band imaging has been obtained, allow- 
ing the selection of line emitting galaxies down to an effective flux 
limit of Fhck 10~^^ erg s~^ cm~^. This has resulted in the largest 
and deepest survey of emission line selected star forming galax- 
ies at z ~ 1, detecting 1517 potential line emitters over an area 
of - 1.4 deg^ in the COSMOS and UDS fields, with 1370 having 
quality multi- wavelength data available in a region of 1.3 deg^. For 
Ha emission line objects this survey probes a co-moving volume 
of ^ 1.8 X 10^ Mpc^ at 2; = 0.84 down to a star formation rate of 
^ 3 M© yr~^ (with an Ah a. — 1 mag extinction correction). 

Photometric redshifts for COSMOS and UKIDSS UDS 
(Mobasher et al. 2007; Cirasuolo et al. 2008) clearly show that the 
majority of the selected emitters are Ha emitters dX z ^ 0.85 with 
a secondary population of [Olll] 5007/H/3 emitters dXz ^ 1.4-1.5. 



Almost 120 emitters were confirmed spectroscopically, from which 
93 are Ha at 2; = 0.84. The contamination within the sample of 
emitters is estimated to be lower than ^ 6%, and the contamina- 
tion within the Ha sample is much lower 0%) based on the cur- 
rent samples. A total of 743 Ha selected emitters (based on their 
photometric and spectroscopic redshifts) was obtained in the two 
fields. These were used to calculate the luminosity function after 
correcting for [Nil] flux contamination, extinction, incompleteness 
and filter profile biases. The morphologies of these emitters were 
also investigated. The main conclusions of this work are: 



• The Ha luminosity function at z 
ted by a Schechter function with 



0.84 found is well fit- 
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-1.92±0.10 



Mpc 



a = -1.65 ± 0.15 and L* = 10^^-^^^° °^ erg s"^ This demon- 
strates a strong evolution in the Ha luminosity function compared 
to lower redshifts and agrees reasonably well with previous smaller 
studies dXz ^ 1. 

• The evolution of the Ha luminosity function can be described 
by an increase in (jf and L*, at least out to 2; ^ 1, with L* then 
continuing to rise up to 2; ~ 2 but (jf peaking around z ^ 1 and 
then decreasing at higher redshifts. 

• The integrated luminosity function is used to estimate the 
cosmic star formation rate density, (psfr) z — 0.84: 0.15 ± 
0.01 M© yr~^ Mpc~^ (corrected for 15% AGN contamination and 
integrated down to 2.5 M© yr~^). 

• An accurate determination of the cosmic evolution of psFR 
has been made using a single star formation tracer (Ha) from 2; = 
to 2; = 2.23. This shows a strong rise up to 2; '-^ 1 followed by a 
flattening out to 2; ^ 2.2. 

• Ha emitters dX z — 0.84 are mostly morphologically classed 
as disks in the rest-frame B-band (82 ± 3%). Irregulars account for 
15 =b 2% of the sample and early-type galaxies are only 3 =b 1%. 
Apparent mergers are a significant fraction of the sample (28ib4%). 

• A strong morphology-Ha luminosity relation is found at 
2;=0.84, with the fraction of irregulars rising steadily with luminos- 
ity and the fraction of disks falling. Mergers/non-mergers present 
the same behaviour, and L* (from the total sample) seems to define 
a critical switch-over luminosity between the two populations. 

• Mergers dominate the bright end of the total Ha luminosity 
function at 2; = 0.84 and ^ 20% of the total psfr is due to their 
activity. 

• A similar morphology-Ha luminosity relation is found at 
lower redshift {z — 0.24), consistent with a simple L* evolution. 

These results suggest that the evolution of both the Ha lumi- 
nosity function and psFR change significantly in nature beyond 
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z ^ 1, and that it is entirely plausible that this is driven by the dif- 
ferent evolutionary behaviour of two different populations of star 
forming galaxies. Out to z ^ 1, the integrated psm at each red- 
shift is produced predominantly by disk galaxies; it is therefore the 
evolution of these disk galaxies, rather than that of major merg- 
ers, which drives the strong decrease in the cosmic star formation 
rate density from redshift one to the current epoch - this is in line 
with other recent results (e.g. Lotz et al. 2008 ). At these redshifts, 
the evolution in psFR arises predominantly from an evolution in 
the characteristic space density of the Ka luminosity function (0* ; 
which evolves by more than an order of magnitude between z — 
and z = 0.84), rather than a strong evolution in L*. With their 
relatively quiescent star formation activity, the disk galaxies dom- 
inate the Ha luminosity function at low luminosities, and are thus 
responsible both for the 0* evolution and for setting the faint-end 
slope of the luminosity function. 

In contrast, irregular and merging galaxies appear to domi- 
nate the Ha luminosity function above L*, at z = 0.84 (and also 
at z = 0.24), being responsible for more than 50% of the non- 
extrapolated psFR (integrated down to 2.5 yr~^). The evolu- 
tion of these systems controls the bright end of luminosity function, 
and thus the cosmic evolution of L*. The continued strong evolu- 
tion of L* between z = 0.84 and z = 2.23, and the decrease in 0*, 
suggests an increasing importance of merger-driven star formation 
activity beyond z ^ 1, especially as the irregular Ha luminosity 
function seems to be very similar to the one found at z — 2.23 by 
Geach et al. (2008) when taking into account an L* evolution. This 
is consistent with recent results such as |Shi et aL] ( |2009| ). The com- 
pleted HiZELS survey will provide statistically significant samples 
of Ha emitters at z = 1.47 and z = 2.23 and so provide a direct 
test of this suggestion. 
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